Reproducing spin lattice models in strongly coupled atom-cavity systems 
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In an array of coupled cavities where the cavities are doped with an atomic V-system, and the 
two excited levels couple to cavity photons of different polarizations, we show how to construct 
various spin models employed in characterizing phenomena in condensed matter physics, such as 
the spin- 1/2 Ising, XX, Heisenberg, and XXZ models. The ability to construct networks of arbitrary 
geometry also allows for the simulation of topological effects. By tuning the number of excitations 
present, the dimension of the spin to be simulated can be controlled, and mixtures of different spin 
types produced. The facility of single-site addressing, the use of only the natural hopping photon 
dynamics without external fields, and the recent experimental advances towards strong coupling, 
makes the prospect of using these arrays as efficient quantum simulators promising. 



Introduction: The burgeoning field of quantum compu- 
tation promises much to the science and technology com- 
munity. While the ability to factor large numbers effi- 
ciently may still be some way off, the advances and poten- 
tial applications brought along with the understanding 
and control of quantum processes, from beautiful manip- 
ulations on minute systems [l| through to coherent many- 
body operations Q, cannot be underestimated. One of 
the first such applications is likely to be the simulation 
of one quantum system with another, more easily manip- 
ulated, quantum system. The most general results have 
been expressed by showing how to simulate one Hamilto- 
nian with another with the help of a series of extremely 
fast single-qubit rotations, breaking the evolution down 
into a sequence of stroboscopic pulses which approximate 
the desired evolution [|[ , which is known as a Trotter de- 
composition. However, in physical systems such as opti- 
cal lattices and ion traps, we possess much more direct 
ways of simulating a variety of different systems, merely 
by adjusting periodic potentials using, for example, glob- 
ally applied lasers [4|, making such simulations feasible 
with current technology. 

Of particular interest are models of the form 



/ ^zZiZj + \ x XiXj + XyYjYj, (1) 

(id) 



where denotes all nearest-neighbour pairs on a lat- 
tice of a particular geometry (typically, a ID chain, or 
2D square lattice) and a is the vector of Pauli matri- 
ces X, Y and Z. There are a number of special cases 
which are commonly examined. For example, the Ising 
model (X z ^ 0) in a transverse magnetic field (B x ^ 0) 
is a simple one-dimensional model which exhibits critical 
properties. Others include the XX (X x = X y and X z =0), 
Heisenberg (A x = X y = X z ) and XXZ (A x = X y ^ X z ). In 
two dimensional lattices, such as the hexagonal lattice, 
simple topological models arise. One possible test-bed for 
these ideas is an optical lattice setup where the natural 
Bose-Hubbard Hamiltonian can be manipulated to pro- 
duce these topological, critical and other effects In 



addition, they are capable of creating three-body terms 
and chiral interactions 

Coupled cavities arrays (CCAs) have been initially 
proposed for the implementation of quantum gates @]. 
Recently, intense interest has arisen from the demon- 
stration that a polaritonic Mott transition and a Bosc- 
Hubbard interaction can be generated in these structures 
0, EH El ■ In the same work it was shown that the Mott 
state could be mapped directly to a spin XX model 
These papers lead to a plethora of studies on various 
properties of CCAs in the direction of many body simu- 
lations [12j , quantum computation [l3| and production of 
photonic entanglement [14]. The study of CCAs provides 
a theoretical framework that can be implemented using a 
variety of technologies such as photonic crystals, toroidal 
microcavities and superconducting qubits [l5|, El, [l7l ]. 
Thus, the aforementioned results are not bound to a spe- 
cific physical system. 

In this paper, the aim is to extend this theoretical 
framework by restricting to the on-resonance, strong cou- 
pling, case and examining how one might enrich the simu- 
lated model by incorporating more complex atomic struc- 
tures within the dopants, and by utilising photons of dif- 
fering polarisations; the goal being to achieve as much 
of the generalised model described in Eqn. (TT]) as possi- 
ble without resorting to a Trotter decomposition, which 
imposes additional experimental difficulties. While such 
decompositions are applicable to the original CCA pro- 
posals [H, [Hi EH , a proposal implementing similar models 
through the rapid switching of a number of off-resonant 
time dependent optical fields followe d up by a Trotter 
expansion has recently been proposed [18;] . Coupled cav- 
ity arrays are capable of single-qubit addressing, so the 
corresponding local magnetic fields B in a spin model 
simulation are readily achieved. The key to creating 
the desired ZZ (which was absent in initial proposals 
@) S3, EJ) interaction is by suitably selecting the de- 
generacy of the energy levels of the dopant atoms. We 
show how an atomic V-system is capable of achieving 
this. This will provide the additional benefit that sim- 
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ply by tuning the number of excitations in the system, 
a large range of different, higher dimensional, spin mod- 
els can be simulated. The possibility of simulating high 
dimensional spins in the presence of strong dissipation 
using constant external fields is also currently being ex- 
amined In the present work, for the case of small 
dissipation, we present a simpler scheme utilizing just the 
natural photon hopping dynamics of CCAs, and no time 
dependent external fields or detunings. 

Atomic V System: We start by considering an array 
of cavities, placed on the vertices of an arbitrary lattice 
(typically, we consider a regular lattice such as a ID chain 
or 2D plane). Each cavity is doped with a single system 
(which we refer to as an atom) , whose energy level struc- 
ture is that of a ground state, \g) and two degenerate 
excited states |A) and \B), depicted in Fig. [1] Within 
each lattice site, the Hamiltonian takes the form 

H int = (aa» + 66+ + \A)(A\ + \B)(B\) 
A A \A)(A\ + A B \B)(B\ 
+g(\A) {g\®a+\g) (A\ ® a*) 
+g {\B) (g\®b + \g) (B\®tf), 

where a+ and create photons of orthogonal polarisa- 
tions, and are those responsible for promoting the ground 
state of the atom to the excited states \A) and \B) re- 
spectively. Henceforth, we assume that the atomic lev- 
els and the cavity are on resonance (i.e. the characteris- 
tic frequency of the cavity is equal to the frequency of 
the atomic transitions of the ground state to the excited 
states; = As = 0). The strength g represents the 
strength of the coupling between the cavity and the atom. 

In the basis \i/j,Na,Nb), we can calculate that the 
(unnormalised) on-site eigenvectors are 



y° Sn ) = ^/s~^i\A,n-l,S-ri)-y/n\B,n,S-n-l) 
■ \pn |A, n — 1, S — n) + V S — n \B, n, S — n — 1) 
±V~S \g, n, S — n) 



S.n 



with energies Su>o and Sujq ± gVS respectively (see Fig. 
[TJ). Na and Nb are the number of a and b photons in 
the cavity, and is the state of the atom. Here, n is an 
integer index (0 to S) which enumerates the basis within 
the manifold containing S excitations. 

Let us assume that we are working at unit filling frac- 
tion, so we expect one excitation per lattice site, meaning 
that only the states 

|0) = (|A,0,0)-| 5 ,1,0»/V2 
|1) = (|B,0,0)-| 5 ,0,1»/V2 

are populated. This arises from the observation that 
there is an energy penalty of U — (2 — \2)g for mov- 
ing from one excitation per lattice site to having two 
excitations in one site, and none in another. 




FIG. 1: The atomic V-system on resonance with a cavity. 
There are two orthogonal photon types, <r and 6^ , which only 
cause transitions to a single level (A or B) from the ground 
state. This gives rise to a non-linear internal structure on 
each site. 

The individual cavities are coupled together by an in- 
teraction 

Hho P = J a (a\a l+1 + a t a\ +1 ) + J b (b\b l+1 + b t b\ +1 ), 

where J a , Jb <C U correspond to the hopping strengths 
for the two different polarizations of photons between 
neighbouring cavities [3]. The effect of the coupling can 
be studied by applying perturbation theory (to the sec- 
ond order) to a pair of neighbouring sites, using the for- 
mula 



H cS = ]T 

a,fc£{0.1} 2 



\b) (a\J2 



(b\ H hop \n)(n\H hop | a) 
E - E„ 



where are all possible eigenvectors involving 2 excita- 
tions on 1 site, and none on the other. Calculating the 
relevant matrix elements in the |0),|1) basis we find the 
effective interaction Hamiltonian 

H eS = -B z {t®Z+Z®t)-\ z Z®Z-\ x {XX+YY) (2) 



where 



k = 



A, 



31 
Wg 

9 
32g 



Jh) 



B, 



5 

9JgJ b 

16g 



{ Jl - Jl) 



and we have ignored the term k1 which simply con- 
tributes a global phase. The local magnetic fields can be 
manipulated by applying local Stark fields of our own, 
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FIG. 2: A comparison of the ground state energy between a 
simulation of the full system (dots) and the prediction from 
perturbation theory (solid line) for 4 cavities, doped with an 
average of 1 excitation per site. The chosen parameters are 
g — 1CP 3 , J a = 10 -5 . A phase transition occurs at J a = Jb 
between the |0)® 4 and |0}® 4 ground states. The energies have 
been scaled to remove the shift of 4u>o — 4<?. To observe other 
phases, such as the one in the XXZ model requires compen- 
sation of the B z term by external fields. 

thereby leaving an XXZ Hamiltonian where the coeffi- 
cients X z and A^ are independently tunable (at manu- 
facture of the device). A comparison of the theoretical 
prediction and an exact diagonalization are depicted in 
Fig. (TJJ) . A degree of tunability of the Hamiltonian can 
be introduced at run-time by varying the detunings of 
the atomic transitions. However, one must remain in the 
regime where the detuning is small so that the pertur- 
bative expansion still holds, which restricts the range of 
variation. 

Generalised Model: Our hopping terms, with strengths 
J a and Jb, effectively describe transmission of photons 
(between cavities) through a birefringent crystal with fast 
and slow axes aligned with the directions a and b. In an 
optical lattice, one can rotate these axes by applying a 
Raman transition to the tunnelling potential. In CCAs, 
the ability to apply this rotation is dependant on the 
particular realisation under consideration. In a setting 
where the cavities are connected by optical fibres, such 
as fibre-coupled micro-toroidal cavities [l5[ , these optical 
fibres represent the birefringent material that we require, 
and the optical axes (c and d) can be aligned indepen- 
dently from the directions defined by the atomic tran- 
sitions (a and b). Moreover, the degree of birefringence 
{J a / Jb) and the orientation can potentially be tuned dur- 
ing the experiment by applying an electric field perpen- 
dicular to the fibre, and making use of the Kerr effect, 
rather than having to initialise all of these properties at 
the point of manufacture. In Circuit QED and photonic 
crystal realisations, however, the hopping comes directly 
from the overlap of the wavefunctions of the individual 
sites [HI, [l7j . which are thus directly connected to the 
a I b basis, and it seems unlikely that these will support 
this generalisation. In cases where this rotation can be 
achieved, the two sets of axes are unitarily related, 



FIG. 3: In a hexagonal lattice, with spin-1/2 located at each 
vertex, and the indicated couplings, topological effects arise. 



and the simulated Hamiltonian is changed to (V 
V)H c g(V <g> VY . While this generates a variety of dif- 
ferent terms, for example XiY 2 +YiX%, we are unable to 
realise the fully anisotropic model XYZ. 

One very useful simulation that is introduced due to 
this rotation is that of the hexagonal lattice ||. At one 
limit, this yields the toric code [20j, and in another region 
yields non-Abelian anyons with the aid of an external 
magnetic field. It is readily formed by setting Jb = 0, 
which implies that X x = 0, and then rotating, along set 
directions, the remaining term ZZ into XX and YY as 
required (see Fig. [3]) . 

Within the optical lattice community, the possibility 
of setting X x ~ X z = has been explored with a view to 
eliminating two-body terms, so the leading order of per- 
turbation theory gives three-body interactions. Armed 
with this toolbox, one couldgenerate many interesting 
effects such as chiral terms [7|. In optical lattices, this 
possibility is achieved by using a Feshbach resonance, 
such that the collisional energies U can be tuned arbitrar- 
ily. In the present system, in order to set X x = X z = 0, 
one requires J a = Jb = i.e. the spins are not coupled, 
and so three body terms cannot arise. We might hope 
to mimic the effect of Feshbach resonances by introduc- 
ing a detuning between the atom and the cavity, which 
would serve to shift the energy levels. However, in or- 
der to maintain the system's integrity, such a detuning 
should be Aa,b <C g, in which case the shift in energies 
is unable to entirely cancel the X z term. 

Higher Spin Models: Unlike the simple two-level 
dopant considered in [9j], changing the average number 
of excitations per site influences the Hamiltonian that 
is simulated. If there is an average of S excitations per 
site, where S is an integer, then there are 5+1 ground 



states, 



= v 



for n = to S, enabling the simulation 
of a spin-iS* particle. Again, there is an energy barrier 

of U = (2VS - VS + T - y/S~^l)g ~ gS^' 2 to having 
any number other than exactly S excitations on each lat- 
tice site, so the ground state is the Mott phase for small 
J/U. All of these models can simulate a Hamiltonian of 
the form in Eqn. @, except with differing coupling co- 
efficients, where the spin operators take on the form of 
the generalised SU(2) X, Y and Z rotations respectively 
for the spin ^S. For example, with 2 excitations per site, 
we realise an array of qutrits interacting through a form 
described by Eqn. where X, Y and Z are replaced 
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by the equivalent qutrit operators, 
1 



Jx = | 1 1 
1 



Jz = —i[Jx, Jy] and 
124\/2 , 2 



Jy 




A, 



7.9 
123 



J Q 2 + J 2 ) S 2 = 



53 



7V2.g V a 



Jl + Jb) A, 



l23V2J a J b 

7g 



Again, further refinements can be incorporated by imple- 
menting the polarization rotations due to the presence of 
a birefringent material. If the rotation is described by 
the 2x2 unitary, 



V 



-ie(n :I X+nyY+n z Z) 



then the effective Hamiltonian is rotated by 



V 



The functional form of the coupling constants for ar- 
bitrary S can be calculated, but is pathological. We 
note, however, that the leading order matrix elements arc 
0{\fS). This significantly adds to the diversity of mod- 
els that can be efficiently simulated in this simple model, 
just by changing the number of excitations present in the 
initial state of the system. 

Non-integer filling: Given that an integer number of 
excitations, S, per lattice site describes spin-iS 1 parti- 
cles, a non-integer value of average excitations per site 
potentially describes a blend of different types of parti- 
cles. Consider the general case where the filling fraction 
isS' + Z, 0</<l. The minimum energy configuration 
is for a mixture of particles of spin and ^(S+l) in the 
ratio (1 — /) : /. The analysis of first-order perturbation 
theory on Hh ap , yields, for the low energy dynamics, a 
swapping of the particles between the sites, governed by 
the effective Hamiltonian 



#eff * S ,i) *S+1,, + 1 



4(5 + 1) 



(j q V(* + i)(j + i) 



+ J b y/(S-i + l)(S-j) 



S+l,i 



S+M+l 



S,j 



S,j+1 



which should be symmetrised for the possibility where 
the higher spin particle starts on the left. As already 
discussed, to first order, there is no interaction between 
particles of the same type. While we are unaware of 
a physical phenomena that this simulates, it completes 
the analysis of the system in the on-resonance case, and 
demonstrates the potential that coupled cavity arrays 
possess. 

Conclusions: We have described a scheme to realize a 
family of spin systems in an array of coupled cavities. By 



introducing a V-configuration to the dopants, the range 
of nearest-neighbour Hamiltonians that can be simulated 
is vastly enhanced. With an integer average of S exci- 
tations per site, we simulate nearest- neighbour spin- iS 
interactions. For S — 1, the spin-i model allows us to 
reproduce the Heisenberg, XX and XXZ models as well 
as those that exhibit both phase transitions and topolog- 
ical features. In the case of a non-integer filling fraction, 
we simulate a mixture of two particle types interacting. 
The resultant strong spin-spin coupling and the individ- 
ual addressability of the separated cavity-atom systems 
make this approach a promising step towards the real- 
ization of quantum simulators for many-body spin prob- 
lems. Since completing this work, we have become aware 
of other work which has considered the same V-system 
[23| |. which just considered the case of S = 1, recovering 
the same results presented here. 

The results presented here for the simulation of spin- 
as' particles are exact, up to terms 0(J 3 /g 2 ) and in the 
absence of decoherence (J = max(J a , Jb)). The primary 
causes of decoherence are photon loss from the cavities 
and spontaneous emission from the atoms, whose rates 
are given by k and 7 respectively. If decoherence is 
present in the system, our results remain valid while the 
corresponding rates are dwarfed by the effective hopping 
rates, which requires \[SJ 2 jg ^> max(-\/SK, 7). 

Note that even though this scheme is not especially 
robust against decoherence, such losses cause the sys- 
tem to leave the computational subspace and are thus 
detectable in the final measurement steps. In order to 
more successfully combat decoherence, one must utilise 
a scheme where the states of interest are ground or dark 
states rather than excited states. Work is progressing 
in that direction, with results to date requiring the use 
of a complex scheme employing constant external fields, 
an elaborate detuning configuration and weakly coupled 
cavities Nevertheless, it may be possible to find in- 
teresting regimes within our model where quantum phe- 
nomena persist, even in the presence of decoherence. For 
example, in [2l| . it is described how entanglement can 
persist in the steady state between a pair of noisy cavities 
when coupled through a third, pumped, cavity. Although 
this work makes no reference to how such a scheme might 
scale, or what information might usefully be extracted, 
it suggests that further investigation is warranted. The 
case of non- integer filling fraction is, in fact, more robust 
to decoherence because it only utilises first-order per- 
turbation theory, and hence we work in a regime where 
g >• JS max(ySK, 7). For the case of circuit QED 
recently g/ max(K, 7) ~ 400 has been reported [161 ]. 

Another intriguing case to study is the atomic V- 
system in the off-resonant case, and see how the be- 
haviour of the two different photon types mimics those 
of two-species or single species spinor Bose Condensates 
(see, for example, [221]). which should be different in na- 
ture to the non-integer fractional filling discussed here (it 
has the potential to allow particles to change type). 
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